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ABSTRACT 

The chemical evolution history of a galaxy hides clues about how it formed and has 
been changing through time. We have studied the chemical evolution history of the 
Milky Way (MW) and Andromeda (M31) to find which are common features in the 
chemical evolution of disc galaxies as well as which are galaxy-dependent. We use a 
semi-analytic multi-zone chemical evolution model. Such models have succeeded in 
explaining the mean trends of the observed chemical properties in these two Local 
Group spiral galaxies with similar mass and morphology. Our results suggest that 
while the evolution of the MW and M31 shares general similarities, differences in the 
formation history are required to explain the observations in detail. In particular, we 
found that the observed higher metallicity in the M31 halo can be explained by either 
a) a higher halo star formation efficiency or b) a larger reservoir of infalling halo gas 
with a longer halo formation phase. These two different pictures would lead to a) a 
higher [0/Fe] at low metallicities or b) younger stellar populations in the M31 halo, 
respectively. Both pictures result in a more massive stellar halo in M31, which suggests 
a possible correlation between the halo metallicity and its stellar mass. 

Key words: galaxies: chemical evolution - galaxies: halos - galaxies: formation - 
galaxies: evolution 



1 INTRODUCTION 

The chemical properties of galaxies hide clues about their 
formatio n and evolution. Semi-analytic chemic al evolution 
models iTalbot fc Arnetd Il97ll: iTinslevI Il98d) have suc- 
ceeded in explaining the mean trends of galactic systems by 
numerically solving a set of equations governing the simpli- 
fied evolution of the chemical elements as they cycle through 
gas and stars. One strength of these models is that they typ- 
ically have the fewest number of free parameters, making 
convergence to a smaller set of solutions more likely. 

Strong constraints can be placed on chemical evo- 
lution models only by contrasting them against a com- 
prehensive set of observed properties. Since the most de- 
tailed observational data are generally available for the 
Milky Way (e.g. iFreeman fc Bland-Hawth orn .2Q02i), suc- 
cessful agreement between model predictions and these ob- 
served properties has been obtained by several studies in 
the past, which help us to understand the formation his- 
tory of the Milky Way (MW). Some models have focused on 
the evolution of the chemical abundances both in the solar 
neighbourhood and in the whole disc, adopting a frame- 
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work in which the MW has been built-up inside-out by 
mean s of a single accretion event (e . g . iMatteucci fc Francoid 
Il989l : IPagel fc Tautvaisiend Il995l) . Others have used an 
early infall of gas to explain the thick disc formation, fol- 
lowed by a slower infall to form the thin disc (e.g. Chi- 
appini, Matteucci fc Gratton 1997). Several studies have 
paid particular attention to the chemical evolution of a 
larger range of elements (e.g. Timmes, Woosley fc Weaver 
1995; .Goswami fc Pr antzos 2000; Alibes, Labay fc Canal 
2001) or focused on the features of the Metallicity Distri- 
bution Function (MDF ) in the solar neighbourhood (e.g. 
iFenner fc GibsonI l2003l) . However, an implicit assumption 
remains unanswered: is the MW a typical spiral? 



Andromeda (M3 1) is the closest spiral to our MW (e.g. 
Ivan den Berghll2003l) . Previous theoretical studies of the 
chemi cal evolution of M31 have been done bv lDiaz fc Tosil 
("1984") and MoUa, Ferrini fc Diaz (1996), both emphasising 
the evolution of the M31 disc. They concluded that M31 has 
a formation history and chemical evolution similar to that 
of the MW. More recently, however, there has been consid- 
erable observational progress in the study of both galaxies. 
A striking difference between M31 and the MW is that the 
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metallicity of the M31 halo (([Fe/H]) ^ -0.5) ^ is sigriif- 
ican tly higher than its M W analogue (([Fe/H]) ^ —1.8, 
e.g. iRvan fc NorrisI Il99ll ). as revealed by many recent 
studies (e.g. Holland, Fahlman & Richer 1996; Dur- 



rell, Harris, fc Pritchet 2001: JSa^edin^^^M^Du^ne 
20011: iFereuson fc JohnsonI I2001I: iReitzel fc Guhathakurta 



iFerguson et al l 20021: IWorthev fc Espafial l2003l:_ 
i ni et al.l l200a iBrown et alJ l2003t iReitzel et alJ 



20021: 

BcUazzh^^^^JT^^, 

2004; Rjc^SllioS . 

It is therefore timely to attempt the construction of a 
chemical evolution model for both the MW and M31, using 
the same framework. Such an attempt may be helpful in 
highlighting the common features in the chemical evolution 
of spiral galaxies (at least in these two spirals) and those 
which remain galaxy- dependent. 

In Section 2 we describe our multi-zone chemical evo- 
lution model, and in Section 3 we present the results of our 
MW model. In Section 4 we show the results for M31. Fi- 
nally, our results are discussed in Section 5. 



2 THE MODEL 

In this study we use a semi-analytic multi-zone chemical 
evolution model for a spiral galaxy. Th is model is b ased 
on GEtool jFenner fc G ibson 2003; Gibso n et ai]|2003ll . We 
follow the chemical evolution of the halo and disc of M31 
and the MW under the assumption that both galaxies are 
formed from two phases of gas infall. The first infall episode 
corresponds to the halo build-up, and the second to the 
inside-out formation of the disc. Similar formalisms have 
been successful in modeling the chemical evolution of the 
solar nei ghbourhood (e.g. C hiappini, Matteucci & Graf- 
ton 1997: Ichang et al.lll999l : Alibes, Labay & Canal 2001). 
We assume that halo stars were born in a burst induced 
by the collapse of a single proto-galactic cloud (Eggen, 
Lynden-Bell & Sandage 1962) or by t he multiple merger of 
building-blocks llSearle fc Zinnlll97d iBekki fc Chibal 120011 
iBrook et all2003ll . The disc is assumed to be built-up inside- 
out by the smooth accretion of gas on a longer timescale. 
Observations of HI High Velocity Clouds (H VCs) that ap- 
pear to be currently falling onto the MW fe.g. lPutman et all 
l2003l and references therein) may provide evidence for such 
gas infall. Recently, HVCs have also been detected in the 
M31 neighbourhood, though their interpretation as infalling 
clouds is debated ( Thil kcr ot al. 2004). 

We monitor the face-on projected properties of the halo 
and disc components. While this geometrical simplification 
is suitable for approximating the flat disc, it is less appro- 
priate for the halo, whose shape is roughly spherical rather 
than disc-like. However, we consider this choice acceptable 
in a simplified model of the chemical evolution of a spiral 
galaxy. We follow the chemical evolution of several inde- 
pendent rings, 2 kpc wide, out to a galactocentric radius 
R — lORd, where Rd is the disc scale-length. We also ignore 
the bulge component, because we are interested in the rela- 
tively outer region (i? > 4 kpc). Each ring is a single zone 
onto which gas falls, without exchange of matter between 
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the rings. We trace the chemical evolution of each zone indi- 
vidually. In our model, we assume that the a ge of the galax y 
is t„ow = 13 Gyr. The basic equations (e.g. lTinslevlll98'oll . 
in a zone at a radius R, for the evolution of the gas surface 
density Ej,i(_R, t) of an element i are written as follows: 

^g.i{R, t) ~ 

- il){R,t)X^{R,t) + I ^{R,t~T^) 



X Yi{m,Zt-r^)'^^^dm + k 



Mp 



Mb„ 



+ (1-fc) / i>{R,t-Trrr)Y,{m,Zt-^r,J 

■4>{R,t-Tm)Y,{m,Zt-,„,) 

MBma^ 



X f^dm-\'Xi^ , (t)h(R)e-h 
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-(t-tdelay> 

+ X,,„^(t)d(i?)e ^.<«) . 



(1) 



max , 

is the 



Here X,{R, t) = is the mass fraction for the element 

i; il'{R, t) is the star formation rate (SFR); tp(m) is the Initial 
Mass Function (IMF) with the mass range Mmin - M, 
Tm is the lifetime of a star with mass m; Yi{m, Zt— 
stellar yield of the element from a star of mass m, age r„i 
and metallicity Zt_r„ . The first term describes the depletion 
of the element i which is locked-up in newly formed stars. 
The second and the fourth terms show the contribution of 
mass loss from low and intermediate mass stars. The third 
term describes the contribution from Type la SNe (SNela). 
The contribution f r om SN ela is calculated as suggested in 
ICreggio fc Renzinl lll983ll . where fc, Mb^^^, MBma^^ fJ^min, 
/i, /(/i), are defined. The fifth term shows the contribu- 
tion from Type II SNe (SNell). The sixth and the seventh 
terms represent the infalling halo and disc gas, respectively. 

The Kroupa, Tout fc Gilmore (1993) IMF is used here. 
We have chosen a lower mass limit of Mmin = 0.08 Mq 
and imposed an upper mass limit of Mmax = 60 Mq in 
order to avoid the overproduction of oxygen and recover 
the observed trend of [O/Fe] at low metallicity at the solar 
neighbourhood in the MW. Such IMF upper mass limit is 
currently loosely constrained by stellar formation and evo- 
lution models. Yet, these stellar models, and the yields they 
provide, are one of the most important features in galactic 
chemical evolution, although questions remain concerning 
the precise composition of stellar ejecta, due to the uncer- 
tain role played by processes including mass loss, rotation, 
fall-back, and the location of the mass cut, which separates 
the remnant fr om the ejected mat e rial in SNe. The SNell 
yields are from lWooslev^fc Weaver! il995l). We have halved 
the iron yields shown in'Wooslcv & WeaveiJ 1^9^, as sug- 
gested by,Tim mcs et al. (1995). The SNela yields are from 
Iwanwt^^^^J (|W99|) . The metallicity-dependent yields of 
Renzini fc Volil (|l981j) have been used for stars in the range 



1-8 Mq. The lifetimes of stars as a function of mas s and 
metallicity have been taken from lSchaller et al.l lll992h . 
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Table 1. The parameters for the MW and M31 models. 





a 




Rd 






^delay 




bd 




i^d 






(M© pc-2) 


(kpc) 


(Mq pc-2) 


(Gyr) 


(Gyr) 


(Gyr) 








MW 


2 


6 


3.0 


48 


0.1 


1.0 


2.0 


1.25 


0.125 


0.03 


M31a 


2 


6 


5.5 


46 


0.1 


1.0 


2.0 


1.25 


0.125 


0.03 


M31b 


2 


6 


5.5 


46 


0.1 


1.0 


0.7 


0.50 


12.5 


0.08 


M31c 


2 


57 


5.5 


46 


0.5 


6.0 


0.1 


0.10 


0.125 


0.08 



2.1 Infall 

The infall rate during the halo and disc phases is simply 
assumed to decline exponentially, as seen by the adopted 
sixth and seventh terms in equation (1). The evolution of 
the total^ surface mass density, Et(ii, t), is described by: 



(2) 



^^5^^^ = h{R)e^ + d(R)e 

Here, the first term describes the infall rate in the halo 
phase. The infall time-scale in the halo phase, r^, is assumed 
to be independent of radius, for simplicity. The infall of disc 
gas starts with a delay of tdeiay, as seen in the second term. 
The time-scale of the disc infall depends on radius as follows: 

R 



Td{R) ^ ad + hd-^^Gyr. 



(3) 



The values for the constants and hd are free parameters. 
The infall coefficients h[R) and d{R) are chosen in order 
to reproduce the present-day total surface density in the 
disc, 'Lt,d{RTtriow), and halo, Y.t,h{R,tnow), respectively, as 
follows (e.g. [Timmes et al. 1995^1 : 

h{R) — I^t,h{R,tnow) 

( ^ ~i 



X<Th 



exp 



(tnow \ 
Th ) 



d{R) — I^t,d{R,tnou 



xlraiR) 



1 — exp 



tnow idelay 

MR) ^ 



(4) 



(5) 



The infalling halo gas has been assumed to be of primor- 
dial composition. On the other hand, it is unlikely that the 
accreting gas has primordial abundance at a later epoch, 
since even low density inter-galactic medium, such as the 
Lyman a forest, has a significant amount of metals (e.g. 
ICowie fc Songalialll998l) . and it is known that the HVCs in 
the MW, which may be infall ing gas clouds, have metallic- 
ities between 0.1 and 0.3 Zq iSembach et alJboOSi) . There- 
fore, we assume that the gas accreting onto the disc is pre- 
enriched. The level of pre-enrichment can be loosely con- 
strained from the observed metallicity of Galactic HVCs. 
We simply assume that the metallicity of the infalling disc 
material is Zi„fauiR,t) = Z{R,t) if Z{R,t) < Zi„faii,max = 
0.3 Zq, where Z{R,t) is the metallicity of the gas at 
the radius 7? and the time t, otherwise Zi„faii[R,t) = 
Zinfaii,max = 0.3 Zq . The abuudauce pattern of the in- 
falling disc gas is further unknown parameter. Following the 



^ Hereafter, by "total" density we mean the sum of the stellar 
and gas densities. 



above simple assumption, we set the infalling disc gas, at a 
given galactocentric radius R and time t, to have the same 
abundance pattern as the ISM at R and t. This guaran- 
tees the smooth evolution of the gas abundance and of the 
abundance patterns in each radial bin. 



2.2 Disc and halo surface density profiles 

We adopt the following exponential profile for the present- 
day total surface density of the disc component: 



^t,d{R,tnow) — '^t,d{RQ,t now J c 



(6) 



Here Rq — 8 kpc, which is the galactocentric distance of the 
Sun within the MW. The same definition of Rq = 8 kpc is 
applied to the M31 models. 

The surface density profiles of the MW and M31 disc 
are different. We have chosen a scale-length of Rd = 3.0 kpc 
for th e MW (e.g. Robin, C reze & Mohan 1992: lRuphv et alJ 
'1996'; Frcudcnrcich '1998^ and Ra = 5.5 kpc for M31 
LWalte r bos fc Kcimicul I .1988.) . For the MW, we assum e 
T,t,d{RQ,t„ou,) = 48 Mq pc"^ jKuiiken fc Gilmord [l99lh . 
We adopt Et,d(J?0, tnom) = 46 Mq pc"^ in M31, such that 
the mass of the M31 disc is simi lar to that of the MW disc 
(f» 10" Mq as inl^eema^dii)- 

We adopt a modified Hubble law for the present-day 
total surface density profile of the halo iBinnev fc Tremaind 
Il987l): 



T,t,h{R, tn 



1 -f {R/Rq 



(7) 



Here we set a = 2. This corresponds to a volume halo density 
profile: 



pt,h{r,t„ow) = 



Pt,ho 



3/2 ' 



(8) 



with pt,ho = '2^^'^Pt,h{Re,tnom) and T,t,h„ = 2/?Qpt,ho- 
For the MW, we assume T,t,h{RQ,tnow) = 6 Mq pc~^ 
which yields a present-day stellar surface density at Rq of 
^*,h{RQ,tnow) = 1.3 Mq pc~^ for our model (see Section 
3), which is consistent with the observed halo stellar density 
at the solar radius p*,;i(7?Q, tnom) = 5.7 x 10* Mq kpc""^ as 
estimated by Preston, Shectman fc Beers (1991). 

The assumption of q — 2, which implies Pi,^h oc r~^, 
agrees with a recent analysis (Zibetti, White fc Brinkmann 
2003) of the halo emission for a sample of ~ 1000 edge-on 
disc galaxies within the Sloan Digital Sky Survey (SDSS). 
This result is similar to the density profile oc r~^'^ which 
has been suggested for the MW stellar halo JChiba fc Beer j 
l2000l:IChiba fc Beerdl200ll: Sakamoto, Chiba fc Beers 2003). 
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2.3 Star Formation Rate 

We assume that the halo star formation (SF) happens on 
a short time-scale because of a rapid infall event associated 
with the collapse of a single massive proto-galactic cloud 

iEgg en et al. 19623 or multiple mergers of building-blocks 
Searle fc Zinnlll97al . The disc SF is assumed to be a more 
quiescent phenomenon, and l ikely to be driven by the spiral 
arms fe.g. lWvse fc Silklll989ll . Therefore, we adopt a differ- 
ent SF law for each component. 

The adopted halo SFR is described as: 

where Uh is the star formation efficiency (SFE) in the halo. 
Therefore, the halo SFR follows a Schmidt law with ex- 
ponent 1.5 (e.g. |Kennicut3[l99i). The adopted halo SFE 
is Vh ~ 0.125, which is approxi matively half of the value 
{uh = 0.25 ± 0.07) suggested by iKennicu tt (1998) and is 
chosen to give the best fit to the observed halo MDF at 
the solar neighbourhood (see Section 3). Stars born before 
tdeiay, wheu the disc phase starts (Section 2.1), are hereafter 
labelled as "halo stars". 

The adopted disc SFR is written as: 

MR,t) = '^dSaiR,tf % Mo Gyr-i pc-^ (10) 

where Vd is the SFE in the disc. This formulation 
JWvse fc Silklll989ll reflects the assumption that SF in the 
disc is triggered by the compression of the ISM by spiral 
arms. The efficiency factor Ud is a free parameter. We have 
found that Vd affects both the present-day gas fraction and 
the disc MDF. The value Ud = 0.03 is used in our MW model 
to reproduce the observed gas density profile of the MW disc 
and the observed MDF at the solar neighbourhood. 




[Fc/H] I (Gyr) 

Figure 1. Comparison between the results of MW model 
(solid lines) and observations of a) the MDF in the MW at 
the solar neighbourhood (closed boxes with error-bars, Ko- 
toneva et al. 2002), b) halo MDF at R = Rq (Ryan & Nor- 
ris 1991, closed boxes with statistical Poissonian error-bars), c) 
[O/Fe] and [Fe/H] for the stars observed at the solar neigh- 
bourhood (Carretta et al. 2000, who included reanalysis of Sne- 
den et al. 1991, Tomkin et al. 1992, Kraft et al. 1992, Edvards- 
son et al. 1993; Gratton et al. 2003; Bensby, Feltzing & Lund- 
strom 2003; Cayrel et al. 2003), and d) the present-day SFR at 
the solar neighbourhood as summarised in Rana (1991, the shaded 
region). 



3 THE MILKY WAY MODEL 

Using the multi-zone model described in the previous sec- 
tion, we construct a model which most closely reproduces 
the known observational properties of the MW. The adopted 
parameters are summarised in Table 1 (MW model). 

Fig. la compares the model MDF at the solar 
neighbourhood with the observed MDF of K dwarfs 
jKotoneva et alJl2002f) . Here we have assumed a K dwarf 
mass range of 0.8 - 1.4 Mq and convolved the MDF of the 
model with a Gaussian error function with a = 0.15 dex, 
consistent with the known empirical uncertainties of the 
observational data. There is good agreement between the 
model and the data. In Fig. lb, the predicted halo MDF 
at the solar radius is shown. Our MDF includes all stars 
still on the main-sequence at the present-time. The MDF 
has been convolve d with a Gaus sian error function with 
a = 0.25 dex (e.g. lAsplundll2003ll . Our model is in agree- 
ment with the observed hal o MDF at the solar neighbour- 
hood jRvan fc Norrislll99J) . The evolution of [O/Fe] as a 
function of [Fe/H] for the gas component in the solar neigh- 
bourhood is shown in Fig. Ic. The model result is consistent 
with the trend of [O/Fe] and [Fe/H] observed in local stars. 
Fig. Id displays the predicted SF history at the solar neigh- 
bourhood. T he SF history has been estimated by severa l 
authors (e.g. iBerteUi fc NasillioOll : iHernandez et alJl200Cl) . 



Unfortunately, such observational estimates are not defined 
well enough to provide useful constraints on a chemical evo- 
lution model. Nevertheless, Fig. Id demonstrates that our 
MW model is consistent wit h the broad r ange of the esti- 
mated SFR as summarised in iRanal 1^^). 

Left panels of Fig. 2 show the predicted radial distri- 
bution of hydrogen column density (upper) and the oxygen 
(middle) and sulfur (bottom) abundances of the gas at the 
present-day, and compare them to the observations. The ob- 
served hydrogen surface density i s obtained by summing the 
surface densities of H2 and HI l|Dameiri993fl '^. The result 
of our MW model is compatible with the observed hydro- 
gen distribution at the inner radii within the uncertainties 
(« 50%) of the observed values However, the hydrogen 
surface density at the outer radii is slightly overestimated 
when compared with the observations. 

The predicted radial abundance profiles of oxygen and 
sulfur also reproduce the abundances observed in HII regions 
JVilchez fc Esteban|[l99S lAffierbach et alJll997li . although 

^ The HI surface d ensity profile has been recently confirmed by 
iNakanishi fc Sofud j£o03). 

■* The contribution of HII should also be considered. However, 
so far no reliable measurement has been achieved, since it is not 
straightforward. 
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Figure 2. Radial distributions of tlie liydrogen surface density 
(upper) and of tlie oxygen (middle) and sulfur (bottom) abun- 
dances at the present-day for the MW model (left) and the 
M31 models (right). In the left panels, the solid lines are for 
the MW model results. In the right panels, solid, dotted and 
dashed lines display the M31a, M31b and MSlc models, respec- 
tively. The observed distributions of the total surface density of 
hydrogen for the MW are from Dame (1993, open boxes), and 
those for M31 are from Dame et al. (1993, open boxes) and 
Loinard et al. (1999, closed boxes). Observations of abundances in 
HII regions in the MW from Vilchez & Esteban (1996) and Afher- 
bach et al. (1997) (closed and open boxes, respectively) are shown 
in the left panels. Those in M31 from Dennefeld & Kunth (1981) 
and Blair et al. (1982) (closed and open boxes, respectively) arc 
shown in the right panels. 



at the outer radii the model slightly overestimates the sulfur 
abundances. 

Our semi-analytic model, whose parameter values sum- 
marised in Table 1, satisfies the general MW observational 
constraints. We now use the same framework to study the 
chemical evolution of M31. 



4 THE M31 MODELS 

The M31 models employ for the disc a different surface den- 
sity profile from that of the MW model. Since the halo den- 
sity profile of M31 is unknown, we adopt the same halo pro- 
file used in the MW model (see Section 2.2). The SF history 
of halo and disc is described by the parameters listed in Ta- 
ble 1. In Section 4.1, we first show a M31 model which adopts 
the same parameter set of our MW model. This model fails 
to reproduce some crucial features observed in M31. There- 
fore, in Sections 4.2 and 4.3, we present two models able to 
better explain the key observations. 



4.1 M31a model: MW analogue 

First, we construct a M31 model (M31a) with the same pa- 
rameter values of our MW model. These parameters are 
summarised in Table 1. The right panels of Fig. 2 show the 
radial distributions of the hydrogen surface density and the 
radial profiles of the oxygen and sulfur abundances of the 
gas phase. The radial profiles of oxygen and sulfur abun- 
dances are reproduced within the observational errors. The 
M31a model results in too high hydrogen surface density 
when compared with the data. This is because the observed 
hydrogen surface density in M31 is smaller than in the MW. 

Fig. 3 shows the radial profile for the mean [Fe/H] 
(([Fe/H])) of main-sequence (MS) stars for the M31a model. 
Hereafter we simply call ([Fe/H]) the "mean metallic- 
ity". Model results a re compared with observations by 
iBellazzini et al.l l|2003)^- We have chosen as reference fields 
those which are mostly disc- (or halo-) dominated, with 
estimated halo- (or disc-) contamination aroun d or lower 
than 10% (see Table 1 in IBellazzini et alJ l200:^: the disc- 
dominated fields (G287, G119, G33, G76, G322 and G272) 
lie at deprojected^ galactocentric distances i? « 8 kpc, 
R X 12 kpc, « 13 kpc, « 14 kpc, R ^ 15 kpc and 
i? ~ 18 kpc, respectively; the halo-dominated fields (G319, 
Gil, G351, G219 and Gl) lie at projected galactocentric 
distances i? ~ 16 kpc, R~ 17 kpc, R~ 19 kpc, ii « 20 kpc 
and i? ~ 34 kpc, respectively. Fig. 3 shows that M31a is 
in broad agreement with the observed mean metallicity in 
the disc-dominated fields, though the metallicity gradient 
is slightly steeper than the observed one (Table 2). On the 
other hand, the mean metallicity of the halo in M31a is too 
low, compared to the data. 

To overcome the failure of M31a, we explored the pa- 
rameter space of the M31 models which could explain the 
observational properties of M31, and found two viable solu- 
tions. In the following, we present these two models. 

4.2 M31b model: our best model 

The model parameters of M31b are summarised in Table 1. 
Compared with the M31a model, i.e. the MW analogue, this 
model has a shorter time-scale for the infall of disc gas (i.e. 
smaller ad and 6^) and higher disc and halo SFE. We found 
that a combination of higher disc SFE and shorter time- 
scale for the disc infall leads to better agreement with the 
observed radial profile of the hydrogen surface density. In 
addition, a higher SFE in the halo phase leads to a more 
metal-rich halo, and we adopt a SFE 100 times higher in 
the halo phase than the M31a model. 

In M31b the radial profile of the oxygen and sulfur 
abundances reproduces the observational data to roughly 
the same degree of the M31a model (Fig. 2). On the other 
hand, the higher disc SFE leads to stronger SF and there- 
fore larger depletion of gas, and lower present-day hydrogen 
surface density than in M31a. Consequently, the result of 

^ The observed MDFs are derived from red giant branch (RGB) 
stars, whereas the model produces the MDFs of MS stars, since it 
is difficult to construct MDFs of RGB stars within our framework. 
We assume that this inconsistency does not affect our comparison 
significantly. 

^ The inclination angle of the M31 disc is ijv/3i ~ 12.5°. 
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Figure 3. Present-day radial profile of the mean stellar metallic- 
ity of both disc (upper panel) and halo (lower panel). Solid, dot- 
ted and dashed lines represent the results of the M31a, M31b and 
M31c models, respectively. The mean metallicities from Bellazz- 
ini et al. (2003) are also shown (closed boxes with 1 a dispersion 
of their MDF). 
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Figure 4. The MDFs of the M31b and M31c models (dotted 
and dashed lines, respectively). The MDFs observed from Bel- 
lazzini ot al. (2003) are also shown (closed boxes). 



the M31b model is in better agreement with the observed 
radial profile of the hydrogen surface density (Fig. 2). The 
hydrogen surface density profiles in both the M31b model 
and the observations pe ak at a radius of ~ 8 kpc with a 
value of ~ 4.5 pc~^ iLoinard et alJlT999^ . However, like 
in the MW model, the hydrogen surface density at the outer 
radii is over estimated in M31b when compared with obser- 
vations from lDame et all lll993l) . We have tried to reproduce 
the observed hydrogen surface density at the outer radii by 
changing various parameters, however, we are unable to find 
a better parameter set than the M31b model. This might 
suggest that, to explain the low gas surface density at the 
outer radii in M31, another mechanism which is not included 
in our semi-analytic model is required. 

The additional benefit of the shorter time-scale of the 
disc infall in M31b is its smaller SFR at the present- 
day. Since M31b has a SFR which peaks at an earlier 
epoch than M31a, M31b predicts a lower present-day SFR 
(« 1.5 M0yr"^) in the disc (4 - 20 kpc) than M31a 
(« 2.4 MQyr~^). The observational estimates of the present- 
day SFR in M31 still have a large uncertainty. For exam- 
ple, WiUiams (2002) estimated the SFR for the last ^ 5 
Gyr b etween roughly 2 and 20 Mq yr~^, while Iwilliamj 
l|2003h inferred that the total SFR for the M31 disc has been 
« 1 M© yr~^ over the past 60 Myr. This value is higher 
than the SFR of ~ 0.2 M(7) yr~^ es timated from Ha and Far 
Infra- Red luminosities (jPevereux e t al. 1994). In addition, 
all these v alues could be affected by systematic errors (e.g. 
Isdil 120031) . Nevertherless, these observations suggest that 
M31 has a lower SFR than the MW, and M31b is consistent 
with this trend. 

Fig. 3 shows improved agreement between the results 
of the M31b model and the observations of the mean metal- 
licity for both disc- and halo-dominated fields. In the disc, 
M31b has a shallower metallicity gradient than M31a, and 
better reproduces the observed gradient (Table 2). This is 
mainly due to the small bd we adopted (Table 1). 

Fig. 4 directly c ompares the MDFs of M31b with 
the MDF obtained in iBellazzini et alJ l^2003^ in different 
fields. The model MDFs are derived from main-sequence 
stars, and convolved with a Gaussian error function with 
(T = 0.25 dex. There is qualitative agreement between the 
MDFs of the M31b model and those observed, especially 
in the disc-dominated fields. The main exception is the Gl 
halo-dominated field, whose observed MDF is more metal- 
poor and bimodal. The MDFs of the G351 and G219 halo- 
dominated fields show narrower peaks than those of the 
model. Differences in the shape of the observed halo MDFs 
among the different fields suggest possible inhomogeneities 
in the M31 halo (see also JBellazzhii ct al. 2003.) . 

In the framework of our semi-analytic models, M31b 
gives the best fit for the observational constraints available 
in M31. In Section 5, we discuss the formation process of 
M31 which is implied by these results. 

4.3 M31c model: another possible model 

The result of the M31b model suggests that a stronger halo 
SF can produce a metal-rich halo. However, applying the 
higher SFE is not the only way to induce stronger SF. This 
can also be achieved by employing a larger reservoir of in- 
falling halo gas, without increasing the SFE. Here we discuss 
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Table 2. Radial gradients for the mean stellar metallicity of the 
MDFs in the M31 models. The results for the MW model and for 
the fields observed in Bellazzini et al. (2003) are also presented 
as reference. 





(A([Fe//f])/AR) « 


(A([Fe//f])/AR) b 


MW 


-0.058 


-0.021 


MSlobs 


-0.026 


-0.008 


M31a 


-0.041 


-0.021 


M31b 


-0.021 


-0.003 


M31c 


-0.014 


-0.009 



This value (dex/kpc) refers to the range 4-16 kpc for the 
MW model and 4-20 kpc for the M31 models. 
^ This value (dex/kpc) refers to the range 4-30 kpc in the 
halo of the MW and M31 models. 

the M31c model, which assumes a much larger present-day 
halo surface density and a longer (6 Gyr) halo phase (Ta- 
ble 1). Again, to explain the observed hydrogen surface den- 
sity, M31c uses a shorter disc infall time-scale, and a weaker 
radial dependence, which also leads to a present-day SFR 
similar to M31b. 

Fig. 2. shows the predicted radial distributions of oxy- 
gen and sulfur abundances which are in agreement with the 
observational data to roughly the same degree of the other 
models. The hydrogen surface densi t y pro file of the M31c 
model is consistent with lDame et all lll993l) in the inner re- 
gion of the disc, R < 12 kpc. However, M31c has a broader 
peak at a larger radius than M31b. Consequently, the model 
overestimates the surface density in the outer region, R> 12 
kpc. Again, this result might suggest that additional physi- 
cal mechanisms operate, to explain the observed low hydro- 
gen surface density at the outer radii. 

Fig. 3 shows that M31c also reproduces the observed 
radial distributions of the metallicity, similarly to M31b. 
The metallicity gradient in M31c is flatter than M31b at 
the inner radii as a consequence of the weaker radial depen- 
dence of the disc infall time-scale (i.e. a smaller bd)- In the 
halo fields, M31c also has systematically higher metallicity 
than M31a and a steeper gradient than M31b. Although the 
gradient is in better agreement with the observations than 
M31b (Table 2), M31c has lower mean metallicities than ob- 
served, especially at the outer radii. In Fig. 4, the predicted 
MDFs at different radii are displayed. The MDFs of M31c 
are also in good agreement with those observed, especially 
in the disc. On the other hand, the halo MDFs show a more 
pronounced metal-poor tail than in M31b. 

Thus, the M31c model also leads to acceptable results, 
although M31b is in better agreement with the observations. 
The next section discusses the formation history of M31 im- 
plied by the M31b and M31c models. 



5 DISCUSSION 

The results in the previous section have shown that the main 
trends of the chemical evolution history of both the MW 
and M31 can be described using the same framework. We 
have also found that to explain the observations in more 
detail, different sets of parameter values for the formation 



history are required for the MW and M31, respectively. In 
this section, we discuss how these different parameters values 
relate to different galaxy formation histories. 

The previous sections have shown that the main dif- 
ferences between the MW and M31 include the hydrogen 
surface density in the disc and, more significantly, the halo 
metallicity. The smaller hydrogen surface density in M31 
can be explained by a combination of the shorter disc infall 
time-scale and the higher disc SFE. 

The most striking difference is that the mean metallicity 
observed in the M31 halo is about ten times higher than that 
observed in the MW halo, though both galaxies have similar 
mass and morphology. Which is a "typical" halo? Metal-rich 
halos s eem more common than metal-poor ones, as pointed 
out bv lZibetti et al.l ||2004|) who analysed a sample of about 
1000 edge-on disc galaxies within the SDSS. lHarris fc Harris! 
(2001) have also shown that the NGC5128 halo MDF closely 
resembles that in the M31 halo. This similarity is perhaps 
suggestive of a common history in the halo formation of both 
galaxies, despite their different Hubble-type. 

A straightforward way to obtain a more metal-rich halo 
would be by requiring homogeneous pre-enrichment of the 
infalling halo gas with a metallicity Z ~ 0.1 Zq. How- 
ever, this metallicity seems too high, if such pre-enrichment 
occurred homogeneously in the whole universe. For exam- 
ple, although this value is close to the mean metallicity of 
Damped Lyman a systems (PL As), many DLAs have m uch 
lower metallicity JPettinil l20o3 iProchaska et alJl2003l) . If 
DLAs are more evolved systems than the infalling halo gas, 
their metallicity should be higher than that induced by pre- 
enrichment. Therefore, a pre-enrichment of ~ 0.1 Zq may 
not be tenable. 

We found that the observed metal-rich halo of M31 
could be explained by two scenarios without any pre- 
enrichment: a) a higher halo SFE; b) a larger reservoir of 
infalling halo gas with a longer halo phase. Scenario a) might 
be explained in a hierarchical clustering regime, where the 
halo is formed by accretion of building-blocks. Theory pre- 
dicts that gas removal by supernova-driven winds should be 
more effective in lower mass systems, leading to a consequent 
supp r ession of the SF R and a low SFE (e.g. iDekel fc Silkl 
119861: lEfstathioullioOOl) . In addition, observations based on 
the SDSS also suggest that the SFE decreases with de- 
creasing stellar mass in l ow mass galaxies with a stellar 
mass M* < 3 X 10^" Mg jKauffmann et alJl20oj) . There- 
fore, higher SFE in the M31 halo phase would be explained 
if the building-blocks of the M31 halo are systematically 
higher mass systems than those of the MW halo (see also 
iHarris fc Harris 2001). This result again supports the notion 
of diversity in the formation of spiral galaxies which are ap- 
parently similar in mass and morphology and belong to the 
same Local Group. Admittedly, our model does not adhere 
to the hierarchical-clustering scenario in a self-consistent 
manner, however, our chemical evolution models can be in- 
terpreted as the "mean" SF and chemical evolution history 
of the stars which end up at each radius. 

As an observational consequence of our models, we 
found that both scenarios a) and b) predict a more mas- 
sive stellar halo in M31, respectively about 6 and 9 times 
more massive than in the MW, whose stellar halo mass is 
~ 10^ Mq. This result suggests that there might be a cor- 
relation between the halo metallicity and its stellar mass. 
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Using dynamical simulations, Bekki, Harris & Harris (2003) 
show that the stellar halo comes from the outer part of the 
progenitor discs when the bulge is formed by a major merger 
of two spiral galaxies. Based on this, they also predict the 
correlation of the metallicity of the stellar halos and the mass 
of the bulges which were formed by major mergers, since 
larger bulges have the larger progenitors, and progenitor spi- 
ral galaxies should follow the observed mass-metallicity re- 
lation. Although they do not mention a correlation between 
the masses of the bulge and halo, it would naturally be ex- 
pected. Thus, a major merger scenario could explain our 
conclusion. Future observational surveys will better quantify 
the correl ation between the hal o metallicity and its stellar 
mass (e.g. iMouhcine eral]|2003ri . 

It is possible to distinguish scenarios a) and b) through 
observation. First, due to a longer halo phase, scenario b) 
predicts intermediate-age population in the M31 halo. This 
picture could explain the recent evidence of inte rmediate-age 
^popul ation by deep imaging of the M31 halo l|Brown et alJ 

It is worth mentioning that the metallicity gradi- 
ent in the stellar halo is sensitive to the assumed stel- 
lar halo density profile especially in scenario b); a steeper 
density profile leads to a stee per metallicity gradient. 
iPritchet fc van den Berghl il994h suggested that the outer 
halo of M31 can be modeled by a power law surface bright- 
ness profile of / oc _R~*, which is much steeper than what 
we assumed (oc R~^). We found that such a steep profile 
rules out scenario b) to reproduce the fiat metallicity gradi- 
ent observed in the M31 halo. However, it is still difficult to 
accurately measure the halo surface brightness of M31 (e.g. 
IZibetti et alJlioO^ . Thus, more observational estimates of 
the M31 halo surface brightness profile would be an impor- 
tant test for this scenario. 

We also found that the higher halo SFE in scenario a) 
leads to about 0.2 dex higher [O/Fe] when [Fe/H]< -1 dex, 
due to intense halo SF. Although measuring [O/Fe] is a 
hard challenge for the current available instruments, this 
task could be accomplished by the next-generation large- 
aperture optical telescope. 



5.1 Prospect 

The framework we have used can explain the main trends 
in the chemical properties of both the MW and M31. How- 
ever, recent observation s of stellar streams bot h in the MW 
(e.g. iHelmi et al.l Il999|: IChiba fc Beerj l200d: llbata et'al 



20021: iBrook et alJ 



1200; 

'2004^ 



" iNavare^'^Helm^^^ 

and M31 Jlbata et al 



reeman 



20041: iMa iewski et al.' ' 

McConnachie et al 200 3: Mcrr ett et al J 12003: Zuc ker et alJ 
20041: iLewis et al.ll2004ll clearly identify inhomogeneities 
in the chemical and dynamical history of both galaxies, 
which could suggest that a significant fraction of the halo 
stars results from later accretion of satellite galaxies. In the 
light of these recent observations, it will be important to 
study both the MW and M31 in more detail by employing 
a self-consistent chemo-dynamical model (e.g. iBrook et alJ 



l2004f) to trace their interrelated chemical and dynamical 
histories. 
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